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We analyze the entanglement properties of the Bogoliubov vacuum, which is obtained as a second 
order approximation to the ground state of an interacting Bose-Einstein condensate. We work 
on one and two dimensional lattices and study the entanglement between two groups of lattice 
sites as a function of the geometry of the configuration and the strength of the interactions. As 
our measure of entanglement we use the logarithmic negativity, supplemented by an algorithmic 
check for bound entanglement where appropriate. The short-range entanglement is found to grow 
approximately linearly with the group sizes and to be favored by strong interactions. Conversely, 
long range entanglement is favored by relatively weak interactions. Working with periodic boundary 
conditions we find some surprising finite size effects for the very long range entanglement. No 
examples of bound entanglement is found. 

I. INTRODUCTION 

In the recent years there has been a considerable interest in studies of entanglement in quantum distributed systems. 
This is a newly developing interdisciplinary field in which quantum information theory meets atomic and molecular 
physics, quantum optics, condensed matter physics, and quantum statistical physics. There are several areas in which 
the role of entanglement and quantum information in distributed systems may be studied. The first motivation has 
come from the studies of quantum macroscopic and mesoscopic phenomena in atomic physics such as Bose-Einstein 
condensation (BEC) U\&\M The 

group of K. Burnett was perhaps the first to study squeezing and entanglement of 
quasi-particle excitations in trapped Bose condensates, as well as characterization of dynamical quantum states of a 
zero temperature BEC 0,0,0- I n fact, in multicomponent BEC the non- linear interactions can lead to a squeezing 
of the collective atomic spin, opening thus a possibility of applications of BEC for precise frequency measurements. 
This line of research was initiated by Cirac and Zoller and their collaborators 0, Q . 

The new impulse to study entanglement in quantum statistical systems has come from the papers of Amico et al. 
fifty and Nielsen ^lj and his collaborators, who have considered scaling properties of (short range) entanglement close 
to a quant um p hase transition. Follo wing up on these studies, various spin chain models were considered by a number 
of authors [T^ [13, [13, [13, EH [13, [13, EM lla H3, Perhaps the most interesting result obtained so far in this 
context concerns the approach of the Garching group, who has not looked only at two parts of the distributed systems 
by tracing out the rest. On the contrary, Cirac and his colleagues considered localized entanglement of the two parts 
by performing optimal local measurement on the rest of the system Xlm . These authors were able to show that there 
exists an entanglement length diverging at quantum critical points |25| . 

Somewhat independently O'Connor and Wootters has introduced general studies of spin chains and rings looking 
for optimal conditions for ent ang lement. The optimization has been relate to certain nearest-neighbor Hamiltonian; 
first in an approximate sense [2|| I27I |2^ | , but later exactly • 

Apart from the spin systems, harmonic chains and rings have been studied |62j |. As we will describe below, powerful 
theoretical tools are available for these systems. An important advantage over spin systems is that entanglement 
can be studied also between two subsystems containing many sites. In Ref. [3(]|, Audenaert et al. used this fact to 
consider several different subsystems on one-dimensional rings where the coupling between neighboring oscillators is 
only through position operators (i.e. can be interpreted as "springs"). A general lesson learnt was that the ground 
state entanglement between two subsystems decreases rapidly with their mutual separation. Increasing the size of the 
two subsystems in general increases their entanglement, but if the "contact" region between the subsystems is kept 
fixed, a finite limiting value is eventually reached. 

The reason why harmonic chains are relatively easy to treat theoretically is that in a grand canonical description, 
the ground state as well as the thermal states belong to the category of Gaussian states. Gaussian states are very 
nice for investigations from a point of view of entanglement properties, since they are completely characterized by 
their first and second order correlations. This simplification as compared to general infinite dimensional systems 
has allowed a lot of results to be derived. First of all it has been shown that if Alice and Bob have one harmonic 
oscillator mode each and share a Gaussian state, the state is entangled if and only if its partial transpose is not 
positive 0,H2- If Alice has one mode and Bob many, the same conclusion holds. If, however, both parties have 
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more than two modes each, they may share an entangled state which has nonetheless a positive par tial transpose. For 
finite dimensionsal systems the existence of such states was demonstrated by the Horodeckis [33], who have shown 
that these states cannot be used for any entanglement distillation procedures of the kind introduced by Bennett et 
al. |34j. The corresponding distillability problem for Gaussian states was solved with identical result: it was shown 
that Gaussian states are distillable if and only if their partial transpose is not positive 35] . A necessary and sufficient 
entanglement criterion for Gaussian states of two parties was found soon after Q]. All these findings allowed for 
remarkable results concerning classification of Gaussian operations 36] , and in particular the proof of the fact that 
Gaussian states cannot be distilled using Gaussian operations. 

As exemplified by this paper, the abovementioned results also have big practical importance, because the Gaussian 
states of photons and atoms are in many cases the ones that are easily accessible experimentally [33, HH |3^, 0, 
EH l42l l43j. In particular, we will here use some of the machinery for Gaussian states to investigate the naturally 
occurring entanglement in harmonic chains of a very specific kind: The ones that appear in the studies of atoms in 
optical lattices. Ultra-cold bosonic atoms in optical lattices can undergo a super-fluid to Mott insulator transition, 
predicted in Ref. Q (see also H3)> an( ^ observed by Greiner et al. 0|- The Mott insulator state with its regular 
filling is considered as an ideal initial state for quantum information processing 0, |^ , and this has made these 
states the subject of intensive investigations in the recent years. We will here instead consider the super-fluid state 
of the lattice bosonic gas in the intermediate regime where interactions are present, but not completely dominant. 
Then fluctuations and excitations are described by the Bogoliubov-de Gennes equations 49], i.e. are formally just a 
system of coupled harmonic oscillators. We therefore find a similar setting as in Ref. [3(| albeit with a different form 
of the coupling between the oscillators. As we consider also two-dimensional lattices, we are faced with an even wider 
choice of subsystems. We have chosen to focus on subsystems that each consist of a string of contiguous sites and to 
vary the size, the separation, and (in the two-dimensional case) the relative orientation of the two strings. Finally, 
we also vary the one remaining physical parameter, namely the ratio between the energy associated with tunneling 
between sites and the mean-field interaction energy. 

Let us briefly summarize our results. As the separation between the two subsystems grows, their entanglement 
decreases and eventually disappears entirely. This is very similar to the results for the "spring" -chains of Ref. [30| . 
Using periodic boundary conditions a peculiar finite size effect occurs and the entanglement does not only depend 
on the distance between the closest parts of the two subsystems, even when this is very well-defined. For the one- 
dimensional strings we consider, the dependence of the entanglement on their (equal) length is essentially linear up 
to moderate separations. This holds even when the relative orientation of the strings is changed. The dependence 
on the interaction strength is such that the entanglement at small separations is increased for increased interactions 
while the entanglement at large separations is decreased. For a given separation, there is thus an optimal ratio of 
tunneling and interaction energy with respect to maximizing the entanglement. Somewhat surprisingly, we find no 
examples of bound entanglement. 

The rest of the paper is organized in the following way. First we give some necessary background of quantum 
information theory in Sec. |HJ In Sec. IIIII we then calculate the relevant correlations of the system. In Sec. IIVI we 
apply the theory to this input and we present and analyze our results. Finally, in Sec. we discuss our findings and 
their relations to similar work. 



II. GAUSSIAN STATES AND ENTANGLEMENT 



In general, the study of entanglement in continuous variable systems is quite demanding. The important subset 
Gaussian states is however much easier to treat and for these states the entanglement properties are fairly well 
understood 0,H(|. The Gaussian states are well known from the field of quantum optics, since important classes like 
thermal states, coherent states, and squeezed states are all Gaussian. In the field of atom optics, Gaussian states also 
appear, albeit slightly less naturally since one tends to prefer a description in terms of Fock states, i.e. states with a 
definite particle number, when dealing with massive particles. This point is discussed further in Sec. [3 below. For 
the moment, however, let us just assume that we have M modes described by "position" and "momentum" operators 
with canonical commutation relations 

[x u ,p v >] = i5 uv > , [x v ,x v i] = [p v ,p v >]=Q. (1) 

To simplify notation, it is useful to refer to x v and p u as Tiv and r 21^+1, respectively. Then Eq. can be written as 

[rv.rv] = iJ vv :, (2) 

where J is the so-called symplectic matrix 

l" 
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Gaussian states can now be denned as states for which the Wigner characteristic function |50( is Gaussian: 

Xw{€) = (exp (i£ T r)) = exp (i$, T d + jr T jr j . (4) 

In Eq. J2J|, d = (r) is the average displacement while 

7„„ = 2R<(r„-d I ,)(»v - d v .)) ( 5 ) 

is the covariance matrix. Note that a Gaussian state is fully determined by specifying <i and 7. In fact, since 
displacements can always be removed by local operations, only 7 is important for the entanglement properties of a 
state and in the following subsections we summarize how information about the entanglement is extracted from it. 



A. Logarithmic negativity 

An important problem in quantum information theory is how to meaningfully quantify the amount of entanglement 
in a state. Viewing entanglement as a resource for performing tasks forbidden by classical physics, one seeks measures 
that quantify to which extend a state is useful for a certain task or how many standard resources it would take to 
create the state. 

Here we apply as our measure of entanglement the logarithmic negativity |5l| . En- While this measure is by no 
means perfect, it does have some very nice properties for our purposes. First of all, the value of En provides an upper 
bound on the efficiency of distillation, i.e., the extraction of maximally entangled states from a larger number of less 
entangled states. Secondly, En is an additive quantity which means that it behaves "naturally" when applied to more 
than one copy of a state. This facilitates comparisons between systems with different numbers of modes. Thirdly, 
and quite importantly, the logarithmic negativity is computable, i.e., there is an efficient way to calculate it for the 
states we are interested in. 

Formally, the logarithmic negativity is defined for any bipartite density operator p as 

^(p)=log 2 ||p TA ||i, (6) 

where p TA is the partially transposed (w.r.t. Alice) density operator and || • ||i denotes the trace norm. For Gaussian 
states one finds |3(| 

2M 

E N =-^min{0,log 2 [\X u (iJj)\] } , (7) 

where Afc(iJ7) denotes the fcth eigenvalue of iJj and 7 is the covariance matrix of /9 Ta . On the level of covariance 
matrices, partial transposition is implemented by reversing all momenta belonging to Alice's subsystem. 



B. Checking for bound entanglement 

One fundamental drawback of En is the fact that a value of zero does not guarantee that the examined state is 
separable: There exist entangled density matrices with positive partial transpose |H,|53- To complete the picture 
and check whether a given state with vanishing En is indeed separable, we apply a qualitative test in the form of a 
non-linear algorithm devised by Giedke et al. |l| . The idea of this test is to start from the covariance matrix of the 
state and then successively create a series of new matrices that all exhibit entanglement if and only if the original state 
does. At each step, simple sufficient (but not necessary) criteria for separability and for entanglement are applied. 
Quite remarkably it can be shown that a definitive answer will come out in a finite number of steps. 



III. FLUCTUATIONS IN THE NON-CONDENSED MODES 



In this section we describe the calculation of the fluctuations in the non-condensed modes for bosonic particles on 
a lattice. We apply the Bogoliubov approximation which will give a good approximation to the ground state as long 
as almost all atoms are condensed. The calculations are rather straightforward, but we give them in some detail for 
completeness. 
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The Bose-Hubbard Hamiltonian 



Our calculations start with the well-known Bosc-Hubbard Hamiltonian 0, [5j| . We will work with d dimensional 
quadratic lattices and have 

£t B h = 2d j a ) a j - J E { a h' + h - c - } + f E a j a h a i ■ ( 8 ) 

J (J J') 3 

The first two terms in Eq. (JBJ describes the possibility for atoms to hop between neighboring sites. They are the 
discrete equivalents of the kinetic energy with 2c? the coordination number (number of nearest neighbors) of the lattice. 
The last term describes the collisional interaction of atoms occupying the same site. The natural modes of the kinetic 
energy are plane waves and we can rewrite 

#BH = J E C l C k £ kin W + ^ff ^ E_ 4, 4 2 % % +k 2 -k 3 ( 9 ) 

k ^1,^2,^3 



with e(fc) = 2J2a=xy (1 — cos[27rfc (T ]) the non-interacting energy associated with wave vector k. N s is the total 
number of sites and the plane wave annihilation operators arc defined via 

^ = 4fE^ 27r " J ^ (io) 



B. The Bogoliubov Approximation 

Starting from Eq. we now apply the Bogoliubov approximation, i.e. we assume that one mode, the condensate, 
is macroscopically populated and develop the Hamiltonian to second order in the fraction of non-condensed particles. 
Formally, we first make the redefinitions 

c - e- l ^ h (c + ^N) (11) 
C £ -» e-'^Cfc for fc^O (12) 

#BH -> -ffBH-ME4 C fe' ( 13 ' ) 
k 

where N is the number of particles and /i = gn with n = N/N s the density of particles. This takes care of the 
macroscopic population of the condensate and its resulting time-dependence, i.e. we now work in a frame where in 
the fully condensed approximation all modes are in their vacuum state and the k = mode defines zero energy. When 
only terms of Eq. © with at least one factor of N are kept, we get the following quadratic Hamiltonian: 

#quad = H +J2 H h ( 14 ) 

fc#0 

where 

Ho = -j9n ^2cjc + cJ4 + c c ) (15) 

H % = (e(fc) + gn) c\c- k + X -gn (c\S_ t + cgc.g) . (16) 

It is well-known how to diagonalize Eq. I|16(l : the normal modes are squeezed combinations of opposite momenta. 
For each pair (fc, —k) we get two bosonic quasi-particle modes, 

c+(fc) = cosh r]£ C£ 4- sinh rjj; c} -. (17) 

c-(fc) = sinhryg ct + cosli77£ c_£, (18) 

with squeezing strength 



^ = J Xe{k l +2 , (19) 
V Ae(Jfc) 
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where A = J/gn quantifies the relative importance of interactions. Energetically, the modes are degenerate and in 
the limit of vanishing interaction they correspond to the usual plane waves. 

We have taken care of all modes except for the condensate itself, the k = mode. If one tries to actually calculate 
the squeezing of the condensate fluctuations it is found to be infinite since e(0) = 0. This expresses the fact that 
the true eigenstates of the system must have a well-defined number of particles and therefore the displaced vacuum 
assumption of Eq. cannot be stationary [54l l55j . The best we can do at this point is to put in by hand a coherent 
state for the condensate, i.e. to assume that the state we study is annihilated by cq (since we have already subtracted 
the macroscopic population). 



C. Fluctuations in a site basis 



We now have a simple description of the system in terms of squeezed momentum eigenstates. Since we are inter- 
ested in the entanglement between different spatial regions of the lattice, we need to calculate the fluctuations and 
correlations of operators describing the atomic field at each lattice site. To use the notation of most of quantum 
information theory literature on Gaussian states, we will use quadrature operators defined as 

a^ + al — at 

For the rex-correlations we find squeezing, 

(xjxj,) - (xj)(xj,) = JL I 1 + 5>os[27rfc • (j-j')] e - 2r >n | , (21) 
and for the pp-correlations we find anti-squeezing 

(PjPyO - (Pj><PjO = I l + E C0S M ■ G-?)]* 2 "* \ • ( 22 ) 



k^O 



In Eqs. I|21I22[I . the "1" is the contribution from the condensate mode, i.e., the part we have put in by hand. There 
are no xp-correlations. 



IV. RESULTS 



In this section, we present our results for the logarithmic negativity of bipartite states. We begin with results for 
ID lattices in Sec. IIV Al then move on to 2D lattices in Sec. IIV Bl 



A. ID lattices 



As we are using periodic boundary conditions, our ID lattice is in fact a ring. We will define the two subsystems 
between which we want to study the entanglement as two sets of contiguous sites, see Fig. ^ Because of the overall 
translational invariance, only the size of the groups, q, and their separation, s, matters. 

In Fig. |21 we plot En as a function of the separation of the two groups for several different group sizes. We note that 
En is generally an increasing function of the group size and a decreasing function of the separation. Small groups 
even become separable already at moderate separations: s > 6 for q = 1 and s > 12 for q = 2. Note that for q > 1 
we cannot conclude this from the vanishing of En alone, but we have to apply the nonlinear algorithm described in 
SecHTBl 

In Fig. the number of sites is N s = 321 and the finiteness of this number naturally is important when the 
separation and/or the group size is comparable to it. In that case, there are two relevant distances from Alice to 
Bob: in the clockwise direction and in the counterclockwise direction. Remarkably, two separations that individually 
(i.e. as the short distance on an infinite ring) would give rise to En = can still result in an entangled state when 
combined. This is evidenced in the insert of Fig. [21 where the vanishing and the revival of En as a function of s for 
q = 3 is plotted on a logarithmic scale. On an infinite ring, we would have En = for all s > 11, but we see that 
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FIG. 1: Definition of the two subsystems. Alice and Bob are both assigned q contiguous sites and the two groups have a 
separation of s sites (between the extreme sites). Since we are using periodic boundary conditions, i.e. working on a ring, there 
are in fact two distances between the groups: s and \N S — 2q — sj + 2. We use the convention that s is the smaller one. 
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FIG. 2: Logarithmic negativity En as a function of group separation s for different group sizes q. In this plot, N 3 = 321 and 
A = J/gn = 20. When both q and s are much smaller than N s , En is a increasing function of q and a decreasing function of 
s. The small separation part (s < 5) can be reasonably well approximated by an exponential with decay constant ~ 0.3. At 
larger separations, the decay is faster. In fact, for small group sizes, En and thus the distillable entanglement drops strictly to 
zero at quite moderate separations: s = 7 for q = 1 and s = 11 for q = 2. The insert shows a logarithmic plot for q = 3. Here 
the negativity vanishes at s = 12, but due to the finiteness of N s , it reappears at s — 76. A similar picture applies for q = 4, 
while for larger qs, finite size effects always keeps En non-zero for N s = 321. 



at s = 76, i.e. at a short separation of 76 and a long of 321 — 2 • 3 — 76 = 138, En again becomes finite. It should of 
course be noted that the value En reaches when the two groups are on opposite sides of the ring is very low (~ 10~ 6 ). 

Apart from the geometrical aspects regarding the lattice and the definition of Alice's and Bob's subsystems, there 
is one underlying physical parameter in the problem, namely A = J/gn. Small values of this parameter lead to a 
higher population of atoms with k ^ or, cquivalently, to more squeezed quasi-particle modes. If we keep n fixed, 
the Bogoliubov approximation will eventually break down and the system will enter a Mott insulator regime |46| . 
Note, however, that since we specify only the ratio J/gn, even values that lead to a high absolute number of excited 
atoms arc not a priori irrelevant since the Bogoliubov approximation still holds for a high enough n. In Fig. [3] we 
compare the curves for q = 3 and 20 < A < 100. At short distances entanglement is clearly favored by a stronger 
non-linearity (low value of A), but perhaps a little surprisingly, beyond separations of about 5 sites, stronger non- 
linearity actually leads to less entanglement. One possible explanation for this phenomenon can be derived from the 
so-called monogamy of entanglement: Alice's and Bob's subsystems are naturally not only entangled with each other 
but also with the remaining sites in the system, in particular with the sites in the gap between them. Since the total 
system is in a pure state, the mixedness of the Alice-Bob system is exactly due to this entanglement. When the short 
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FIG. 3: Dependence of entanglement upon the value of A. The five curves show En as a function of s for q = 3 and five 
different values of A. At short distances, a higher ratio of hopping to non-linearity leads to less entanglement, while at longer 
distances the opposite picture applies. 
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FIG. 4: Assignments of sites to Alice and Bob on the 2D lattice. Note that we again use periodic boundary conditions, i.e. 
although the lattice is drawn as flat, we work on a torus. The "group size" denoted by q gives the number of sites assigned to 
both Alice and Bob. The "center-of-mass separation" denoted by s cm gives the distance between central sites in the two group. 
Finally, the "orientation" denoted by O labels different arrangements of the sites within each group. 

distance entanglement then grows due to a stronger non-linearity, the mixedness of a highly delocalized Alice-Bob 
system must be expected to increase, leaving less room for Alice-Bob entanglement [5(|. To test the plausibility of 
this explanation, we have examined the purity Trp 2 corresponding to Fig. [3] As expected, the purity has a stronger 
dependence on A at large separations than at short separations, but the effect is not very pronounced and we cannot 
rule out entirely different explanations of Fig. [3J 

B. 2D lattices 

On 2D lattices, we have a wide choice of interesting assignments of sites to Alice and Bob. We have chosen to 
focus on the ones shown in Fig. 0] and will investigate to which extend the structure within the groups influences 
the entanglement. To this end, it is natural to first calculate the entanglement between two single sites. In Fig. [5] 
we therefore show the logarithmic negativity as a function of the separation of the sites. We note that to a good 
approximation En is an isotropic function of the separation and that it decays quickly with increasing distance 
between the two sites. Compared to the ID results with the same J/gn ratio, En is now a factor of approximately 
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FIG. 5: Entanglement between two single sites as a function of their separation. The full lattice is 241 sites wide in both 
directions and A = J/gn = 100. 

10 smaller. 

In Fig. we plot our results for the three different "orientations" of the two groups. Instead of plotting directly 
the logarithmic negativity, we first divide En by the number of sites in each group, 

E N {q, s cm ) = ~E N (q, s cm ), (23) 

q 

and then subtract the single site result (Fig. O evaluated at the "center-of-mass distance" (cf. Fig.0J, 

AE N (q, s cm ) = E N (q, s cm ) - E N (1, s cm ). (24) 

This is a sensible rescaling, at least at short distances where AEN(q, s cm ) turns out to be an order of magnitude 
smaller than E N (1, s cm ). It also almost collapses the curves for group sizes q = 3 and q = 5 onto each other: for the 
orientation O=0 to a very high degree, while for 0=1 and 2 some residual entanglement remains. That the orientation 
0=1 where Bob's string is placed as a continuation of Alice's string leads to the highest entanglement at a given s cm 
is not surprising since this is the orientation with the shortest minimal distance between sites in the two groups. 

In ID we found a remarkable inversion in the dependence of En on A as the separation between groups were 
increased. A similar effect is present in 2D as can be seen from Fig. [7\ There we plot En as a function of A for a 
number of different group sizes and separations. All the curves show a maximum in the plotted A range, i.e. in all 
cases the entanglement is optimized at some finite ratio of hopping to non-linearity. As expected from the successful 
rescaling in Fig. the optimal A has only a very weak dependence on the group size. Like in ID, entanglement 
between groups at large distances has a larger optimal A than entanglement at short distances. 

V. DISCUSSION 

First of all, we should comment on the symmetry breaking Bogoliubov approach that we have applied, i.e. on the 
procedure of putting in "by hand" a coherent state for the condensate mode. It is well known that this approach leads 
to identical predictions for excitation frequencies etc. as the symmetry preserving approach that does not assume 
superpositions of different total numbers of particles |5f|. A simple connection between the two can be made by 
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FIG. 6: Results for a 2D lattice with A = 100 and a size of 241 sites in each dimension. Group sizes are q = 3 (open symbols) 
and q = 5 (filled symbols) and the three pairs of curves correspond to the three orientations described in Fig. 2] The rescaled 
quantity AEn is an order of magnitude smaller than the single site logarithmic negativity at separations s cm < 11 indicating 
that the logarithmic negativity is approximately proportional to q and that s cm is a good measure of the distance between the 
groups. 
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FIG. 7: Dependence of En upon A for the orientation O=0 on a 2D lattice (241 sites in each dimension). The dotted curves 
are single site results, q = 1, the full curves are for groups of q = 3 sites each, and the dashed curves are for q — 5. Results for 
4 different separations are shown: s cm = 5, 7, 9 and 11. The three uppermost curves are for s cm = 5, the next three for s cm — 7 
and so on. In the displayed A range, all the curves show a maximum, i.e. En is not a monotone function of A. 



averaging over the phase of the coherent state leaving one with a Poissonian mixture of number states. Instead, our 
results should be seen as describing the entanglement existing between Alice and Bob when they share knowledge 
of the overall phase of the condensate. Of course, this phase cannot be re gard ed as a completely classical piece of 
information as it can only be defined relative to some reference condensate 57]. For recent investigations regarding 
the operational implications of super selection rules see e.g. Ref. |58| and references therein. 

Turning now to our results, the success of the rescaling (Eq. J22J)) displayed in Fig. HJland Fig.0is quite remarkable. 
For the configurations we have considered, the dominant behavior of the logarithmic negativity can be understood 
in terms the single site result evaluated at s cm : it should simply be scaled with the number of sites in each group. 
This scaling is exactly what one would expect from the additivity of En if the sites in the two groups paired up in 
a natural way, but such a pairing is far from obvious when looking at Fig. It is an interesting question for future 
studies to understand more precisely why this scaling applies. One hint can maybe be drawn from recent studies of 
the asymptotic behavior of "spring" -chains [59|. 
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Another interesting conclusion we can draw is that a high degree of squeezing (strong non-linearity) in the system 
does not necessarily increase the entanglement: On the contrary, it will eventually tend to decrease long distance 
entanglement. Similar results have recently been found in related few site systems, both for the ground state jH3 
and in the entanglement dynamics |6l) . It is not clear whether a useful picture of the long distance behavior can be 
derived from the short distance one via the principle of monogamy of entanglement, but it would be very interesting 
to at least study the bounds that exist. Since the Bogoliubov model is both relevant and accessible, we expect our 
results to be useful in such studies, but perhaps the prospects for gaining insight in this way are even greater in 
quantum critical systems [ToL ITU . 

As a perspective, experimental measurement of the entanglement or even implementation of quantum information 
processing in the system would of course be interesting. However, in the formulation given here, access to the full 
covariance matrix of the Alice-Bob system is assumed and a reformulation in terms of more easily accessible quantities 
would probably be necessary. 
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